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Based on the Dirac equations in the two-dimensional n— flux model, we study the interaction 
effects both in nontrivial gapped and gapless Dirac equations with numerical exact diagonalization 
method. In the presence of the nearest and next nearest neighbor interactions: for nontrivial 
gapped Dirac equation, the topological phase is robust and persists in a finite region of the phase 
diagram; while for gapless Dirac equation, charge-density-wave and stripe phases are identified and 
the phase diagram in (Vi, V2) plane is obtained. When the next-next-nearest neighbor interaction is 
further included to gapless Dirac equation, the topological phase expected in the mean-field theory 
is absent. Our results are related to the possibility of dynamically generating topological phase from 
the electronic correlations. 
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I. INTRODUCTION 

Topological insulator (TI) is a new class of time- 
reversal invariant quantum phase, which is topologically 
distinct from the trivial band insulator. The experimen- 
tal discovery of such phase has simulated the study of 
topological properties. It has been shown that TI has 
many exotic physical properties, which make it have po- 
tential applications in the fields of spintronics, quantum 
computing, etcP^. 

Recently due to the rich physics and fundamental 
interest, the correlation effect in the context of TI has 
been explored^. One aspect of such study is to consider 
the interplay between the interaction and the intrinsic 
spin-orbit coupling. Based on the Kane-Mele-Hubbard 
modeP^l many numerical and analytical works are car- 
ried outPES. In these studies consistent results are ob- 
tained, such as: the stability of TI to small interaction; 
the magnetic phase at strong correlation; the existence 
of quantum spin liquid phase; et al. Another aspect of 
studying the correlation effect is to dynamically generate 
spin-orbit coupling from the electronic correlations, re- 
sulting topological Mott insulator (TMIpH. At the mean- 
held level, TMI has been found in systems with different 
geometries22H2l] Since the present studies on TMI are 
limited to mean-field approximation, it is necessary to 
provide direct evidences from exact numerical methods. 

In this paper, based on the Dirac equations in the 
two-dimensional n— flux model, we study the interac- 
tion effects both in nontrivial gapped and gapless Dirac 
equations. With numerical exact diagonalization (ED) 
method, we calculate the Chern number, the fidelity met- 
ric, the energy spectrum and the electrons' distribution 
to characterize different phases. For nontrivial gapped 
Dirac equation, the topological phase persists in a finite 
region of the (Vi, V2) phase diagram. For gapless Dirac 



equation and in the presence of nearest neighbor (NN) 
and next-nearest neighbor (NNN) interactions, charge- 
density-wave (CDW) and stripe phases are identified and 
the phase diagram in (Vi, V2) plane is obtained. However 
when the next-next-nearest neighbor (NNNN) interac- 
tion is further included, the topological phase expected 
in the mean-field theory isn't identified. Our results are 
related to the possibility of dynamically generating topo- 
logical phase from the electronic correlations. 



II. THE MODEL AND METHOD 

We consider a 7r— flux model on square lattice with 
a tight-binding Hamiltoniai i 29 * 30 !, 



Ho = J2u j e tx ^c]c i 



(1) 



J 



where c\ and q are the annihilation and creation opera- 
tors at site iv When the sites % and j are NN neighbors, 
Uj = h and Xi,i+± = 0,Xi,i+v = Ki x . While for the 
NNN case, ty = t 2 and Xi,i+x+y_ = Xi+x,i+ v = 7r«x- The 
phase factors are shown in Figflja), in which a two-site 
unit cell (A, C) can be chosen. In the reciprocal space, 
the Hamiltonian is written as Hq — ^ k -0j c i?o(k)V'k with 
ipk = {ca 1 c c ) t and 

-ffo(k) = 2t\ cos k x cr x — 2t\ cos k y a y — Atz sin k x sink y a y , 

with fJ XtV ^ z the Pauli matrices. The energy spectrum is 
given by 



Ek = ±\J M\ (cos 2 k x + cos 2 k y ) + 1Qt\ sin k x sin k y , 

which is symmetric around zero and has a gap 4|t2| a t the 
two inequivalent Dirac points K = (ir/2,±ir/2). In the 




FIG. 1: (Color online)(a) The hopping amplitudes of Hamil- 
tonian Eq.llJ: nearest neighbor solid links £1, vertical dashed 
links — ti, diagonal links iti in the direction of the arrow, 
(b) The Dirac point in the spectrum when £2 = 0. (c) the 
two degenerate configurations with CDW order and the four 
configurations with the stripe order. 



Chern number of the system is calculated and equals 
1(— 1) for £2 > 0(< 0), which is consistent with 
the analytical result. Alternatively since the system 
possesses inversion symmetry, the topological property 
can be determined by the product of the parites of 
the wave functions at the four time-reversal invariant 
momenta Tj = (riiTr/2)x + (mj7r/2)j/ with n^mi = 
0, ipMS if we select site A of the unit cell as 
the center of inversion, the parity operator writes as 
V[i/) A (r),i[>B(rf,i[>c(r),'il>D(r)] = [rp A (-r) , ipB (-r - 2x - 
2y),ipc(- r ~ 2a;),'0£)(-r — 2y)]. In the momentum space 
it becomes V k = diag(l, e - 2 ^ kl+k2 \ e~ 2ikl , e _2i * 2 ). We 
calculate the product of the parities of the occupied states 
and it is —1 for the topological phased 

In the following of the paper, we use ED to study 
the effect of interactions in Eq.dTI). In this method, the 
wave functions and eigenvalues of several lowest states 
can be directly obtained. To characterize the topologi- 
cal property of the system, we calculate the Chern num- 
ber and fidelity metric of the ground state. In an inter- 
acting system, the Chern number can be defined using 
the twisted boundary phase 61,63, i.e., replacing ki,k2 
with 9\ , 02 in Eq. ([3]) . The fidelity metric g is defined 
as g(V,SV) = I (sv)l V) with the fidelity F(V,SV) = 
\(^o{V)\^o{y+^V) the overlap of the ground-state wave 
functions at V and V + SV^. In the following we set 
ti = 1 as the energy scale and all ED calculations are 
performed on 4 x 4 system. 



vicinity of the two Dirac points, the low-energy Hamilto- 
nian is described by the Dirac equations: 



iOQ 



-2t\k x a x + a2tik v cr z — a4t2<T y 



(2) 



with the valley index a = ±1. In order to include the 
staggered CDW and stripe orders, a four-site unit cell 
(A, B, C, D) can be chosen. Then in the reciprocal space, 
with the reduced Brillouin zone {k : |fea,|,|fey| < 7r/2}, 
the Hamiltonian writes as H' Q — X)k V'k^oWV'k with 
V>k = (ca,c b ,c c ,c d ) t and 



H' (\t) = 2ti cos k x a x ® I + 2ii cos k y a y (g) a y 
—4*2 sin k x sink y a z <8> <J V . 

For *2 7^ the system is topological with gapless 
states associated with the edges traversing the gap. The 
topological phase can be characterized by the Chern 
number, which is defined as^, 



C = 



2-ni 



d 2 kF 12 {k). 



(3) 



BZ 



Here the field strength F 12 (k) = d kl A 2 (k) - d^A^) 
with the Berry connection A 1 ( 2 )(k) = (n(k)|9fe 1(2) |n(k)) 
and |«(k)) the normalized wave function. An analyti- 
cal calculation from Eq.M yields |sgn(i 2 /*i) for each 



\31 



Dirac point, thus the Chern number equals sgn(i 2 /*] 
Using an numerical approach suggested in Ref.(32), the 



III. EFFECT OF INTERACTION TO 
NONTRIVIAL GAPPED DIRAC EQUATION 

We study the effect of interaction in Hamiltonian 
Eq.|l]) and add NN and NNN interactions, 

H^t = } J VinjUj + ^ V 2 n i n j 
(ij) <W» 

By calculating the Chern number of the ground state 
of the total Hamiltonian Hq + Hi nt , the phase diagram 
is obtained. As shown in Fig(2j there are three different 
regions: /: the trivial insulator with CDW order; 77: the 
topological phase with Chern number 1; III: the phase 
with four-fold degenerate ground state. 

At Vi = V 2 = 0, the system is in a topological phase 
with quantized Hall conductance a xy — Ce 2 /H with the 
Chern number \C\ = 1. As V\ is increased, the system 
is driven to a trivial insulator with CDW order. The 
Chern number change from 1 to at the critical value 
Uf = 0.84. In the energy spectrum as shown in FigMa), 
for Vi < Vi the energy of the ground state is gapped 
from other states. For Vi > Vf the gap reopens, but the 
ground state becomes nearly two-fold degenerate. Fur- 
ther calculations of the distribution of the particles show 
that they correspond to the two different configurations 
of staggered CDW order. In FigMb), the three lowest 
eigencnergies of the system is plotted near the transi- 
tional point. The states denoted by the red dashed lines 




FIG. 2: (Color online)The phase diagram in (VijVb) plane. 
Here £2 = 0.1 at which the non-interacting system is topo- 
logical. Region I: trivial insulator with CDW order. Region 
II: topological phase with Chern number 1. Region III: phase 
with four-fold degenerate ground state. The size of Region II 
shrinks as £2 decreases. 
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FIG. 3: (Color online)The eigenenergies of several lowest 
states and the fidelity metric, (a), (b): Vi = 0; (c), (d): 
Vi = 0.5; (e), (f), (g): Vi = 3. (b), (d), (f) and (g) are the 
zooming views near the transition points of the corresponding 
left figures, which are marked by the circles in (a), (c) and 
(e). Here £ 2 = 0.1. In (c), (d), (e), (f) and (g), at small V 2 
the red dashed lines correspond to the two states with CDW 
order (if not split it is two- fold degenerate) . 



are with CDW order and they are split due to the pres- 
ence of the hopping terms. The figure clearly show the 
crossing between the states, which means the disappear- 
ance of the topological phase. The obtained critical value 
is consistent with that from the Chern number. To char- 
acterize the topological phase transition, we also calcu- 
late the fidelity metric g and the curve shows a jump at 
the critical interaction. 

In Fig|3j we also show the energy spectrum and the 
fidelity metric g at V\ = 0.5,3, which are representative 
cuts in the phase diagram Fig[2j In all cuts the results are 
consistent with those from calculating the Chern number. 
At V\ = 0.5, the system at V% = is topological. The 
NNN interaction drives the topological system to a phase 
with four-fold degeneracy. At V\ = 3 the system at V 2 = 
is a trivial insulator with CDW order. Surprisingly the 
inclusion of NNN interaction can recover the topological 
property in a finite region. 

To understand the above results, we present a qual- 
itative explanation from a mean-field viewpoint. The in- 
teraction can be decoupled in the direct and exchange 
channel, 



Uirij « rn{rij) + {ni)rij - {ni){rij) 



C i C j\ C j C i) C i C j 



If the CDW and stripe orders are considered, we have the 
ansatz («,-) = 1/2 + p(-l) l * +i y + ^(-l) 1 *. Let the value 
of (c[cj) is St for NN and A for NNN, then this procedure 
yields a mean-field Hamiltonian which writes, 

H' (k) = 2(ti - StVi) cos k x a x ® I 

+ 2(£i — StVi) cos kyOy <8> a y 

— 4(i2 — XVs) sin k x sin k y a z <g) a y 

~ Ap(Vx - V 2 )a z <g> I + AvV 2 (t z ®a z + C 

with the constant C* = N{Vl + v ^ + 2N[(Vi - V 2 )p 2 + 
V 2 v 2 \ + 2NV 1 5t 2 + 2NV 2 X 2 . Compared to H {k), the 
hopping amplitudes t\ and t 2 are modified and two ad- 
ditional terms related to the CDW and stripe orders ap- 
pear. Then the phase diagram can be qualitatively ex- 
plained by including the two additional terms to i/ (k). 
When only the term p'a z ® I is included, a topological 
phase transition is driven by this term, as shown in Fig. 4 
(a). However when the term v' a z ® a z is also included, 
the topological phase can be recovered and further be 
broken by larger v' . Since generally the value of p'{v') 
increases with V\(V 2 ), the phase diagram is qualitatively 
explained. Although in a more rigorous way the mean- 
field parameters should be determined self-consistently, 
the underlying physics is the same. 



IV. EFFECT OF INTERACTION TO GAPLESS 
DIRAC EQUATION 

Due to the exchange decoupling channel of NNN in- 
teraction in the mean-field approximation, it is suggested 
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FIG. 4: (Color online) (a) The Chern number vs p' at v' — 0. 
(b) The Chern number vs v' at p' = 0.5. 
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FIG. 5: (Color online) (a) The eigenenergies of the six lowest 
states of Eq.Q at Vi = 2. (b), (c) and (d): the zooming views 
near the transition points, which are marked by the circles in 
(a). The red dashed lines are for the two CDW states and the 
black solid lines are for the four stripe states (if splitted each 
is two- fold degenerate). 
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FIG. 6: (Color online)The order parameters p, v vs V2 at 
Vi = 2. The critical points precisely match those in Fig[5] 



that a Dirac semi-metal (SM) can be driven into the 
quantum anomalous Hall (QAH) phase in the presence 
of longer-ranged interaction, realizing TMI. Previous re- 
sults are obtained using mean-field theory, so it is neces- 
sary to study this problem using exact numerical meth- 
ods to give further confirmation. In the following using 
ED method, we study whether the topological phase can 
be generated by interactions in a system with gapless 
Dirac point. 

We drop the t 2 term and firstly consider the effect of 
NN and NNN interactions. The total Hamiltonian under 
study is, 



H 



Tol 



(ij) 



(ij) 



^ V 2 n inj . (4) 

((ij)) 



We calculate the Chern number in the (Vi, V 2 ) plane and 
the values are all zero, showing that no topological phase 
is induced when only the NN and NNN interactions are 
included. 

Then we calculate the six lowest states and study the 
possible phases in the system. Figure [5] shows a represen- 
tative energy spectrum at fixed Vi = 2. For small V2 the 
ground-state is nearly two-fold degenerate and gapped 
from other states, denoting the states as ip\ and tp 2 ■ Two 
new states can be constructed using a linear combina- 
tion: Vi = ^fOi + V2), "02 = 75 ("01 - ^2 ), which cor- 
respond to the two configurations with staggered CDW 
order (FigJTI). The other four eigenstates have two close 
eigenvalues and each is two-fold degenerate. Also by a 
proper linear combinations, four new states can be con- 
structed and each corresponds to one of the four config- 
urations with the stripe order. As V% increases the gap 
decreases and vanishes at a critical value V 2 cl > after which 
the ground state becomes the one with the stripe order 
[see FigJ5lb) for details]. Then the states with CDW or- 
der have higher energies. From another critical value V£ 2 
the split between the two CDW states begins to decrease 
and vanishes soon [see Figl5tc)]. As V2 further increases, 
the split between the four stripe states begins to decrease 
at the critical value V 2 c3 . Meanwhile the gap between the 
CDW and stripe states vanishes. After that the gap re- 
opens and there is no split any more. The ground-state 
is exactly four-fold degenerate and with the stripe order. 

For the six lowest states, two of them are with CDW 
order and four are with stripe order. The average density 
on each site can be expressed as (nj) = 1/2 + p(— \y*+ % v 
(CDW order) or (n,) = 1/2 + j/(— 1)**(*») (stripe order). 
In FigJ6] we plot p, v as a function of V2 and find the 
discontinuities at the critical points V 2 c2 , V 2 c3 , which pre- 
cisely match the ones in the energy spectrum. At V 2 c2 , 
p jumps to zero while at V 2 c3 , v jumps to nearly 1/2. 
We scan in the (Vi,V 2 ) plane and give the critical lines 
determined by V 2 cl , V 2 c2 , V 2 c3 , which is shown in FiglTj 

Then we can conclude the phase diagram when the 
NN and NNN interactions are included in the gapless 
Dirac equation. Basically the ground state is topological 
trivial with CDW or stripe order, which is separated by 



the solid line in FigjTJ In the upper region of the long- 
dashed line, the ground state is with the stripe order and 
are four-fold degenerate. Between the solid and long- 
dashed lines, the ground state are also with the stripe 
order, however the eigen-energies of the four stripe states 
are split to two values, each of which is two-fold degen- 
erate. 




FIG. 7: The phase diagram in (VijVa) plane obtained from 
ED calculations. The black solid line is determined by V£ l - 
The short-dashed line is determined by Vf 2 ■ The long-dashed 
line is determined by V 2 c3 . 

Till now we obtain the phase diagram in the presence 
of NN and NNN interactions and no topological phase 
is found. Next we add NNNN interaction Hnnnn — 
J2(((ij)\) V^UiUj to Eq.(|4l. We calculate the Chern num- 
ber in (Vi, V2) plane at (Afferent fixed V3 and don't find 
any region showing topological phase. 



calcualtions, we obtain the phase diagram in the (V\, V?) 
plane and find the topological phase persist in a finite 
region of the phase diagram. The results show that the 
topological phase is robust to small interactions. To un- 
derstand the phase diagram, we also present a qualitative 
explanation using the mean-field approximation. 

Next we study the possibility of dynamically gener- 
ate topological phase in gapless Dirac equation. When 
only Vi and V2 are present, the system may be in the 
CDW or stripe phase, which are all topological triv- 
ial. When the NNNN interaction is further included, the 
topological phase is still absent in our calculations. 

Compared to the mean-field approximation^, our 
numerical results are different in the following two as- 
pects: 1, when only the NN and NNN interactions are 
present the SM phase is not identified; 2, after the NNNN 
interaction is further included, the QAH effect is absent. 
The deviations may be attributed to the small size on 
which the ED simulations are performed. However it 
should be noted that previous and the present works have 
shown that even for small sizes the topological property 
already can manifest itself welP^l. So our results can 
demonstrate that at least on small sizes the topological 
phase does not be generated following the mechanism 
in the mean-field scheme. Considering the strong quan- 
tum fluctuations in two dimensions which may make the 
mean-field approximation break down, large-scale numer- 
ical simulations are expected to give further clarification 
on this important problem. 
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two-dimensional Dirac equations of the n— flux model. 
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ical phase. Using ED method we calculate the Chern 
number, the energy spectrum, the fidelity metric and the 
electrons' distribution of the ground-state. From these 
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